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Abstract 

From the perspective of perturbation theory, we propose a systematic procedure for the evalu- 
ation of the derivative discontinuity (DD) of the exchange-correlation energy functional in Kohn- 
Sham density functional theory (KS-DFT), wherein the exact DD can in principle be obtained 
by summing up all the perturbation corrections to infinite order. Truncation of the perturbation 
series at low order yields an efficient scheme for obtaining the approximate DD. While the zeroth- 
order theory yields a vanishing DD, the first-order correction to the DD can be expressed as an 
explicit universal functional of the ground-state density and the KS lowest unoccupied molecular 
orbital density, allowing the direct evaluation of the DD in the standard KS method without extra 
computational cost. The fundamental gap can be predicted by adding the estimated DD to the 
KS gap. This scheme is shown to be accurate in the prediction of the fundamental gaps for a wide 
variety of atoms and molecules. 
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Over the past two decades, Kohn-Sham density functional theory (KS-DFT) jlj] has be- 
come one of the most powerful theoretical methods for studying the ground-state properties 
of electronic systems. As the exact exchange-correlation (XC) energy functional E^dp] in 
KS-DFT remains unknown, functionals based on the local density approximation (LDA) 
and generalized gradient approximations (GGAs), have been widely used for large systems, 
due to their computational efficiency and reasonable accuracy. However, owing to their 
qualitative failures in a number of situations {2 5|, resolving these failures at a reasonable 
computational cost continues being the subject of intense research interest. 

The prediction of the fundamental gap Eg has been an important and challenging subject 



in KS-DFT j642l|. For a system of electrons (A^ is an integer) in the presence of an 
external potential Vextij)^ Eg is defined as 



Eg = \{N)-A{N), 



(1) 



where I(A^) = E{N - 1) - E{N) is the vertical ionization potential and A(A^) = E{N) - 
E{N -|- 1) is the vertical electron affinity, with E{N) being the ground-state energy of the 
A^-electron system. Therefore, Eg can be extracted from three KS calculations for the 
ground-state energies of the A^- and (A^±l)-electron systems. However, such multiple energy- 
difference calculations are inapplicable for the prediction of fundamental band gaps of solid- 
state systems ll8|. 

By contrast, the KS gap Aks is defined as the energy difference between the highest 



occupied molecular orbita. 
of the A^-electron system 



22 



OMO) and the lowest unoccupied molecular orbital (LUMO) 



A 



KS 



(2) 



where ej(A^) is the z-th KS orbital energy of the A^-electron system. Therefore, Aks can be 
obtained from only one KS calculation for the KS orbital energies of the A^-electron system. 
Note that Eg is not simply Aks, but is given by 



E, 



Aks + A, 



(3) 



where 



.^c - hm < ^ , . 



5E^c[p] 



N+ri 



5p{r) 



N-ri 



(4) 
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is the derivative discontinuity (DD) of Exc[p\ |23l-l32l|. As the KS gap (even with the exact 



functional) severely underestimates the fundamental gap |2J, |25|, l3l|, the evaluation of the 
DD is tremendously important. Recently, the importance of the DD in the excited-state 33 1 



and time-dependent 



34 



35| properties has also been highlighted. Although several schemes 



have been proposed for calculating the DD, they can be computationally very demanding 
or large systems, due to the use of the nonlocal energy-dependent self-energy operators 



10 



11 



18| or Hartree-Fock operator 



2l| 



In this Letter, we provide a systematic procedure for the evaluation of the DD, based 
on perturbation theory [36j]. The lowest-order estimate of the DD can be expressed as an 
explicit universal (i.e., system-independent) functional of the ground-state density and the 
KS LUMO density, allowing very efficient and accurate calculations of the DD and, via Eq. 
([3]), the fundamental gap in the standard KS method. 
For the exact KS-DFT, I(A^) 



-tN{N) 



22 



23 



28 



37 



and therefore A(A^) = I(A^ + 



1) = —^N+i{,N + 1). Consequently, Eg (see Eq. ([T])) can be expressed as 



(5) 



which is simply the energy difference between the HOMOs of the A^- and (A^-l-l)-electron 
systems U- By subtracting Eq. (EJ from (EJ, A., (see Eq. ©) can be expressed as 



A,e = ejv+i(A^ + l)-eiv+i(A^) 

= (V'Af+l I -^KS I i'N+l) - I -f^KS I i^N+l)- 



(6) 



Here iJxs = {-^^"^ + ^ext(r) + / f^\dr' + ^xc([p]; r)} and V'i(r) are, respectively, the 
KS Hamiltonian and the z-th KS orbital of the A^-electron system, with fa;c([p]; i") being the 



XC potential and p(r) = Yl!i=i being the ground-state density. i/KS = {- 



2m, 



■V2 + 



Vextij) + / T^3^(ir' -|- fxc([p]; I")} and ipiiv) are, respectively, the KS Hamiltonian and the 



2-th KS orbital of the (A^-l-l)-electron system, with Vxc{[p\]'^) being the XC potential and 
p(r) = Yld=i lV'j(r)P being the ground-state density. 

Aiming to compute A^-c (and hence Eg) using only one KS calculation for the A^-electron 
system (e.g., for being apphcable to sohds), we express eAr4.i(A^-|-l) in terms of {ei{N),ipi{j)}, 
based on perturbation theory js^. 

We choose i^KS as the unperturbed Hamiltonian, and suppose that the unperturbed 
energy levels are nondegenerate. Let A be a dimensionless parameter, ranging continuously 



from (no perturbation) to 1 (the full perturbation). Consider the perturbed Hamiltonian 
Hx given by 

Hx = Hks + XHl (7) 

where the perturbation H'^ = / j^^'^X-^dr' + Vxc{[px]',t^) ~ J r ' | ~ '^a;c([p]; i") 
involves px{r) = X^ila^ (filling the orbitals in order of increasing energy). Here 

{'ipii'r)} and {e^(A^ + 1)} are, respectively, the eigenstates and eigenvalues of Hx: 

Hxi'^ir) = efiN + l)i,^ir). (8) 

Eq. (IH]) at A = 1 is simply the KS equation for the (A^+l)-electron system, as it can be 
verified that (r)} and {ei{N + 1)} are, respectively, the eigenstates and eigenvalues of 
Hx=i. Therefore, e^+iiN + 1) = e^=\(^ + !)• 

Writing H'^, V'i^(i"), and e^(A^ + 1) as power series in A, we have 

if; = + Ai7'(i) + A2if'(2) + . . . , (9) 

^P^iT) = V) + AV^fV) + A'^f (r) + ■ ■ • , (10) 
6f(iV + l) = ef) + Ae« + AMn---. (11) 
Inserting Eqs. ([7D, Q, ([lO]), and ([II]) into Eq. gives 

(ifKs + Aif'(°) + X'H'^'^ + . . ■ )(^f ^ + X4'^ + AVf ^ + ■ ■ ■ ) 



(ef + AeSn A^n ■ ■ ■ )(^f + A^« + A^^f ) + . . . ). 



(12) 



Expanding Eq. fl22|) and comparing coefficients of each power of A yields an infinite series of 
simultaneous equations. 

To zeroth order (A°) in Eq. (!22l) . the equation is 

HKs4'\r)=ef4'\r), (13) 

which is simply the KS equation for the A^-electron system (i.e., the unperturbed system). 
We then have ipf\r) = ipi{r) and ef^ = ei{N). Therefore, eN+i{N + 1) = e^^^{N + 1) ^ 
eS+i = eAr+i(A^). Correspondingly, A^c = eN+i{N+l)-eN+i{N) ^ eAr+i(A^)-eAr+i(A^) = 0, 
and Eg = Aks + A^c ~ Aks- Therefore, to obtain a nonvanishing A^c, it is necessary to go 
beyond the zeroth-order theory. 
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To first order (A^) in Eq. ( 122|) (see Supplemental Material '40j]), the first-order correction 
to the orbital energy is 



and the first-order correction to the orbital is 



/,/,(0) I tt'{Q) I ,/.(0)^ 



Note that p,=o(r) = Ef=t' l^^°(r)P = Ef=t' \i^f\^)? = Ef=t' l^.(r)P = p(r) + PL(r), 
where PL(r) = |'?/'Af+i(i")P is the KS LUMO density of the iV-electron system. Consequently, 
we have 

K=o = e'J I /'^r^i ^^' + ^-([P + Pl]; r) - ^;.e([p]; r). (16) 
As eN+i{N + 1) = e^=\(Ar + 1) ^ e^^l, + eii^^, we have 

A,, = e;v+i(A^ + l)-e7v+i(A^) 



r — r 



and Eg = Aks + ^xc ~ Aks + eir+i- Eq- (|33|) is a key result, showing that the DD can be 
approximately expressed as an explicit universal functional of p(r) and PL(r), and can be 
calculated in the standard KS method without extra computational cost. Note that it can 
also be derived from Eq. (Q by assuming ("frozen orbital approximation (FOA)") that 

^ = l,2,3,■■■. (18) 



To second order (A^) in Eq. f l22|) (see Supplemental Material 40|), the second-order 
correction to the orbital energy is 



where 



^'(1) 



dX 



A=0 



+ 



Spxir') J dX 



dr 



A=0 



X 



I ^ ~ 

N+l 



fxc{[p + pL];r,r') 



(20) 



2j]K(^*(r')V^f)(r')) Ur'. 



Here /a;c([p]; r, r') = 5t>j:c([p]; r)/^p(r') is the XC kernel, the asterisk (*) denotes complex 
conjugate, and 3ft(- ■ ■ ) denotes the real part of (■■■)• From Eq. fl5B|) . we have 



-N+l 



E 



(21) 



Correspondingly, e7v+i(A^ + 1) = + 1) e 

e,v+i(A^+l)-ejv+i(iV) 



+ ^Svli + ^Tvli- This gives A^^ 



.(1) 



ejjl^i and Eg 



N+l 

Aks + A 



Aks + (e^'^ 



TV+l ' '^N+l) 



Extending the process further, the HOMO energy of the (A^+l)-electron system can be 
obtained by summing up all the perturbation corrections to infinite order, i.e., eAr+i(A^ + 
1) = e^_^i{N + 1) = X]^^o^7v+i- Therefore, we can, in principle, obtain the exact A^c = 



eN+i{N + 1) - eAr+i(A^) = J2'^=i ^iv+i and the exact Eg = Aks + A^c = Aks + Zlr 



An) 
= 1 ^Af+l- 



For any finite order truncation of the above perturbation series, if two or more unper- 



turbed states share the same energy, degenerate perturbation theory may be needed |36 |. 
Since the concept of the perturbation to the unperturbed Hamiltonian i^KS remains valid, 
this scheme could be extended to estimate the A^c (and hence the Eg) for the degenerate 
cases based on the corresponding degenerate perturbation theory. 

As mentioned previously, the DD needs to be summed to the KS gap to give the funda- 
mental gap. While the DD given by Eq. (j4]) should be the same as that given by Eq. 
for the exact functional, this property may no longer hold true for an approximate func- 
tional^ For example, for a LDA or GGA, while the DD given by Eq. (HI) is shown to vanish 

3, 



21 



301], we emphasize that the DD can be favorably restored by Eq. ([6]), and subse- 



quently approximated by Eq. ( l33l) . Although a more accurate approximation for the DD 
could be pursued by higher-order perturbation theory, we adopt the DD given by Eq. (!33ll 
(i.e., first-order correction) for simplicity. Accordingly, the fundamental gap is predicted by 
summing Eqs. ([2]) and (l33l) in our Aks + ^xc scheme. 
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Here we examine the performance of various schemes in the prediction of the fundamental 



gaps for the FG115 database 



20[ |. which consists of 115 accurate reference values for the 



fundamental gaps of 18 atoms and 97 molecules at their experimental geometries. The 
fundamental gaps are calculated by our Aks + A^^p scheme, the Akr scheme (by Eq. ([2])), 



and the Eg scheme (by Eq. ([S])), using the LDA 4l| and LB94 42(] functionals, and the 
6-311++G(3df,3pd) basis set, with a development version of Q-Chem 3.2 43|. The error for 
each entry is defined as (error = theoretical value — reference value). The notation used 
for characterizing statistical errors is as follows: mean signed errors (MSEs), mean absolute 
errors (MAEs), and root-mean-square (rms) errors. Note that for the Aks ot Aks + ^xc 
scheme, only one KS calculation for the A^-electron system is required (i.e., applicable to 
solids), while for the Eg scheme, which is the Aks + ^xc scheme with A^c being exactly 
calculated by Eq. ([6]) (with no further approximations), two KS calculations for the A^- and 
(A^+l)-electron systems are required (i.e., inapplicable to solids). 

The calculated gaps are plotted against the reference values in Fig. [1] (for LDA) and Fig. [2] 
(for LB94). For both functionals, as the Aks gaps are shown to be vanishingly small (some 
of them are even negative) for the small-gap (smaller than 10 eV) systems, the DDs are 
essential for the accurate prediction of the fundamental gaps. In fact, even for the large-gap 
(larger than 15 eV) systems, the DDs remain significant fractions of the fundamental gaps. 
As shown in Table HI the MAE associated with the Aks + A^c^ Eg scheme is more than 



40|. 



three times smaller than that associated with the Aks scheme 

Due to the use of "FOA" (see Eq. f fTSj) ) in the evaluation of the DD, the Aks + gaps 
tend to be larger than the Eg gaps. For rare gas atoms (e.g.. He, Ne, and Ar), where this 



40|. 



approximation becomes excellent, the Aks + '^xc gaps are very close to the Eg gaps 

For LB94, the Eg gaps are in excellent agreement with the reference values, due to the 
correct asymptote of the LB94 potential, which is a key factor for the accurate prediction of 



the HOMO energies 



28|, 



37 



38 



42| and, via Eq. ([5]), the fundamental gaps. By contrast, for 



LDA, the Eg gaps tend to underestimate the reference values, due to the imbalanced self- 



interaction errors (as the LDA potential is asymptotical 



28 



y incorrect) in the predicted HOMO 



37 



38 



42| . As the Aks + '^xc gaps 



energies of the N- and (A^-l-l)-electron systems 
tend to overestimate the Eg gaps, it appears that there is a fortuitous cancellation of errors 
in the predicted Aks + ^xc gaps, when compared with the reference values. 



In conclusion, we have provided a systematic procedure for the direct evaluation of the 



DD, based on perturbation theory. The lowest-order estimate of the DD is an exphcit uni- 
versal functional of the ground-state density and the KS LUMO density (see Eq. fl33|) ). 
presenting a simple, efficient, and nonempirical scheme for the direct evaluation of the DD 
in the standard KS method. The fundamental gap can be accurately predicted by the sum 
of the KS gap and the estimated DD. The validity and accuracy of this scheme has been 
demonstrated for a wide variety of atoms and molecules, extending the applicability of KS- 
DFT to an area long believed to be beyond its reach. To further improve the accuracy of 
this scheme, a more accurate functional and a more accurate approximation for the DD 
(based on higher-order perturbation theory) should be adopted, although this will neces- 
sarily be somewhat more expensive. Since the concepts of the DD and the perturbation to 
the unperturbed Hamiltonian iJxs are still valid for solid-state systems, this scheme could 
be extended to estimate the DD (a correction to the KS band gap) for solids, where the 
prediction of accurate fundamental band gaps is very challenging for KS-DFT. Work in this 
direction is in progress. 
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SUPPLEMENTAL MATERIAL TO: RESTORATION OF THE DERIVATIVE 
DISCONTINUITY IN KOHN-SHAM DENSITY FUNCTIONAL THEORY: AN EF- 
FICIENT SCHEME FOR ENERGY GAP CORRECTION 



DERIVATIVE DISCONTINUITY FROM THE PERSPECTIVE OF 
PERTURBATION THEORY 

As discussed in our paper, we aim to express eAr+i(A^ + 1) in terms of {ei{N) , ipi{r)} , 



based on perturbation theory 



Here, we provide the full derivation (the notations used 



here are the same as those used in the paper). As Eq. (12) of our paper gives 
(ffKS + Aff'<»> + A^ff'l") + . . . + Arf.f' + + ■■■) 

= (4°' + A.;" + A^<f' + . . . )(*!°' + Av-i" + AVf ' + ••■), 

expanding Eq. ([22]) and comparing coefficients of each power of A yields an infinite series of 
simultaneous equations. 

Zeroth-Order Theory 

To zeroth order (A°) in Eq. fl22l) . the equation is 

^Ks^f (r)=efVf (r), (23) 

which is simply the KS equation for the A^-electron system (i.e., the unperturbed system). 
We then have tpf\r) = ^^(r) and = ei{N). Therefore, eN+i{N + 1) = e^=\(iV + 1) ^ 
e!J+i = ejv+i(iV). Correspondingly, A^^ = ejv+i(A^+l)-ejv+i(A^) ^ ejv+i(A^)-eA,+i(iV) = 0, 
and Eg = Aks + A^c ~ Aks- Therefore, to obtain a nonvanishing A^c, it is necessary to go 
beyond the zeroth-order theory. 

First-Order Theory 

To first order (A^) in Eq. (122|) . the equation is 

H^s4'\r) + V) = V) + V)- (24) 

11 



Taking the inner product of Eq. fl2^ with iljf''^ , 

As i^KS is Hermitian, the first term on the left-hand side is the same as that on the right- 
hand side. Moreover, we have {ipf^ \ V'l"'') = (V'j I V'i) = so the first-order correction to 
the orbital energy is 

e« = (V^f ) I I V^f )) = I H[^, I V^.). (26) 
Rewriting Eq. (12^ . we have 

(i^KS - ef )V^«(r) = - (r). (27) 

Since {ipf\^)} constitute a complete set, il)f'\T) can be expressed as a linear combination 
of them: 

V^f^W = E«f^fW- (28) 

Note that it is unnecessary to include j = i in the sum, for if %1)^^\y) satisfies Eq. (1271) . so 
too does {il)f\v) + ail)f\Y))^ for any constant a, and we can use this freedom to subtract 
off the il)f\r) term so that ^/'j^(r) is normalized at first order in A. Inserting Eq. (I28l) into 
Eq. (l27|l . we have 

E(4°^ - = -(H'^'' - V). (29) 

Taking the inner product of Eq. (l29l) with ijj'"^'' , we get 

(30) 

If = z, the left-hand side is zero, we recover Eq. ( 126|) : If 7^ we get (e[.°^ — e-°^)a^*^ = 
— (^/'^°^ I I ^/'j-^''), or a^*'' = ^^-^ 1^ — (1)^' ^ , so the first-order correction to the orbital is 

/,/,(0) I TJ'{0) I ,/.(0)^ 

^■"V) - E !o _ ' '/'f'(^)- (31) 
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Note that p,=o(r) = l^.'=°(r)P = Ef=t' = EIV M^)? = p(r) + PL(r), 

where PL(r) = iV^Ar+ili")!^ is the KS LUMO density of the A^-electron system. Consequently, 
we have 

^A=o = e'l ^^^Z^^^' + ^-([/^ + -^l]; r) - vUlp]; r). (32) 
As e7v+i(A^ + 1) = e^^%{N + 1) ^ e^^i + ej^^^, we have 

Axe = ejv+i(A^+ 1) -ejv+i(A^) 

/■/ PlWp^,^,^, (33) 



e 



r — r 



+ j PL(r)|wxc([p + PL];r) 



and Eg = Aks + ^xc ~ Aks + Civ+i- Eq- (|33|) is a key result, showing that the DD can be 
approximately expressed as an explicit universal functional of p(r) and PL(r), and can be 
calculated in the standard KS method without extra computational cost. 



Second-Order Theory 

To second order (A^) in Eq. ( I22l) . the equation is 

HKsif'(r) + ff'<»>«)<"(r) + ff'<'>'*f'(r) 

= <!°VfV) + <!'V."V) + f''/'r(r). 

Taking the inner product of Eq. fl34l) with ipf*^ , we get 



(34) 



(35) 



As -f^KS is Hermitian, the first term on the left-hand side is the same as that on the right- 
hand side. Moreover, we have {ipf^ I '^f'^) = 1 {tpf'' \ ipf^^) = (see Eq. (l3Ti) ). so the 
second-order correction to the orbital energy is 

^ Q(Ar)-e,(Ar) + " ^ " 
13 



where 



^'(1) _ 



dX 



A=0 



X 



I r — r 
I r — r 

N+l 



-7T + 



dX 



5pA(r') 
^ + /a;c([p + PL];r,r')| 



A=0 



(37) 



i=l 



Here 



a;cv Lr J i ; 



r, r 



5^^xc([p]; r)/5p(r') is the XC kernel, the asterisk (*) denotes complex 



conjugate, and 9fJ(- ■ ■ ) denotes the real part of 



From Eq. fl36|) . we have 



7W\ 7W\ ^ WN+i ' Wn+i)- 



(38) 



Correspondingly, eAr+i(A^ + 1) = e7v+i(^ + 1) ~ e^^_^_^ + ej^^^ + e^^_^^. This gives 
e;v+i(A^ + 1) - eiv+i(A^) ~ eSili + ej^+i and = Aks + A,, ^ Aks + [t%^ + eg^+J. 
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FIG. 1. Calculated versus reference fundamental gaps for the FG115 database [20|]. The funda- 
mental gaps are calculated by three schemes (see text for details) using the LDA functional. 
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FIG. 2. Same as Fig. [H but using the LB94 functional. 
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TABLES 



TABLE L Statistical errors (in eV) of the 115 fundamental gaps of the FG115 database 20(], 
calculated by three schemes (see text for details) using the LDA and LB94 functionals. 



Error 




LDA 






LB94 




Aks 


Aks + Aj^c 


Eg 


Aks 


Aks + A^c 




MSE 


-7.22 


0.78 


-2.33 


-7.20 


2.73 


0.03 


MAE 


7.22 


2.11 


2.33 


7.20 


2.74 


0.47 


rms 


7.44 


2.45 


2.57 


7.48 


3.13 


0.63 



LDA 



41(1 and LB94 



TABLE IL Fundamental gaps (in eV) of the FG115 database [20], calculated by three schemes (see 
text for details) using the 



42l | functionals. 









LDA 






LB94 




Molecule 


Reference 


Aks 


Aks + A^, 


Eg 


Aks 


Aks + A^, 


Eg 


H (Hydrogen atom) 


12.86 


4.71 


11.05 


9.19 


9.48 


14.26 


11.16 


He (Helium atom) 


27.23 


18.01 


22.01 


21.81 


20.71 


27.60 


27.40 


Li (Lithium atom) 


4.22 


1.13 


5.49 


3.98 


-2.55 


6.13 


4.72 


Be (Beryllium atom) 


9.66 


3.50 


10.44 


7.51 


3.34 


12.05 


8.93 


B (Boron atom) 


7.99 


0.19 


10.30 


6.19 


-0.16 


11.09 


7.74 


C (Carbon atom) 


9.97 


0.08 


13.38 


8.08 


-0.52 


14.28 


9.59 


N (Nitrogen atom) 


14.74 


4.02 


14.02 


11.94 


3.28 


17.01 


13.15 


(Oxygen atom) 


12.14 


0.29 


14.44 


10.27 


-0.28 


17.79 


12.20 


F (Fluorine atom) 


13.98 


0.12 


18.95 


12.13 


-0.81 


22.36 


13.99 


Ne (Neon atom) 


26.91 


17.25 


23.74 


23.40 


18.77 


27.98 


27.69 


Na (Sodium atom) 


4.14 


0.91 


5.44 


3.88 


-1.92 


6.86 


4.93 


Mg (Magnesium atom) 


7.76 


3.39 


8.15 


6.20 


3.68 


9.71 


8.04 


Al (Aluminum atom) 


5.53 


0.09 


6.74 


4.54 


-0.41 


7.26 


5.68 


Si (Silicon atom) 


6.73 


-0.02 


8.37 


5.76 


-0.52 


8.91 


6.52 
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P (Phosphorus atom) 


9.78 


2.16 


9.56 


8.13 


0.65 


10.87 


8.49 


S (Suliur atom) 


8.23 


0.04 


9.85 


7.25 


-0.64 


11.79 


8.34 


CI (Chlorine atom) 


9.30 


-0.07 


12.47 


8.38 


-0.76 


13.29 


9.35 


Ar (Argon atom) 


18.65 


11.23 


15.78 


15.52 


11.85 


18.38 


18.23 


CHg (Methyl radical) 


9.86 


1.79 


8.25 


8.04 


0.67 


11.21 


9.35 


CH4 (Methane) 


15.06 


9.21 


12.19 


11.47 


10.08 


15.55 


14.78 


NH (Imidogen) (""E ) 


13.17 


3.10 


11.72 


11.04 


1.99 


15.71 


12.11 


TVTTT /A • !• 1\ 

NH2 (Ammo radical) 


11.34 


1.97 


10.56 


9.87 


1.17 


14.46 


11.27 


NH3 (Ammonia) 


11.54 


5.64 


9.42 


8.47 


6.59 


12.96 


11.91 


OH (Hydroxyl radical) 


11.27 


0.16 


11.52 


9.79 


-0.63 


16.50 


11.38 


H2O (Water) 


13.35 


6.57 


11.03 


9.85 


7.65 


15.08 


13.66 


HF (Hydrogen fluoride) 


16.91 


9.00 


14.23 


12.50 


10.40 


19.39 


17.10 


S1H3 (Silyl) 


7.95 


1.32 


8.48 


6.80 


0.15 


9.91 


7.49 


SiH4 (Silane) 


14.03 


8.14 


11.36 


10.55 


8.43 


13.69 


13.11 


PH3 (Phosphme) 


11.82 


6.13 


9.63 


8.72 


6.32 


14.08 


11.34 


TT /tt 1 1^ 1 \ 

HjS (Hydrogen sumde) 


11.00 


5.57 


9.69 


8.44 


5.56 


14.12 


11.06 


HCl (Hydrogen sulfide) 


13.36 


7.02 


13.00 


10.49 


7.06 


15.89 


13.30 


C2H2 (Acetylene) 


13.43 


6.78 


14.94 


9.23 


6.70 


16.23 


13.61 


C2H4 (Ethylene) 


12.57 


5.60 


14.41 


9.01 


5.46 


15.42 


12.30 


C2Hg (Ethane) 


13.41 


7.80 


10.79 


10.06 


8.61 


13.93 


13.24 


HCN (Hydrogen cyanide) 


14.31 


7.87 


17.57 


11.53 


7.71 


18.65 


15.17 


CO (Carbon monoxide) 


15.57 


6.87 


18.28 


12.93 


6.54 


19.07 


15.35 


HCO (Formyl radical) 


9.56 


1.31 


12.24 


7.99 


0.97 


12.83 


9.56 


H2CO (Formaldehyde) 


11.56 


3.41 


14.51 


8.94 


3.27 


15.37 


11.86 


CH3OH (Methyl alcohol) 


11.67 


5.75 


9.76 


8.28 


6.70 


13.87 


12.07 


N2 (Nitrogen diatomic) 


17.88 


8.21 


20.51 


15.84 


7.85 


20.94 


17.30 


N2H4 (Hydrazine) 


10.29 


4.43 


8.49 


7.44 


5.33 


11.80 


10.74 


NO (Nitric oxide) 


10.11 


-0.04 


13.79 


8.76 


-0.24 


14.38 


10.15 


O2 (Oxygen diatomic) (^Sg) 


12.52 


2.06 


13.55 


11.42 


1.81 


16.91 


12.44 


H2O2 (Hydrogen peroxide) 


12.65 


4.66 


17.70 


8.99 


4.57 


19.47 


13.08 
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1 2 (r luorme diatomic) 


15.53 


3.42 


20.95 


14.02 


3.24 


21, 


.51 


15.40 


CO2 (Carbon dioxide) 


14.58 


8.41 


13.44 


11.90 


8.51 


20, 


.97 


15.64 


P2 (Pnospnorus diatomic) 


10.19 


3.71 


11.25 


9.37 


3.46 


11, 


.58 


9.96 


S2 (Sulfur diatomic) ( Sg) 


7.96 


1.05 


9.64 


7.40 


0.55 


9.1 


65 


7.75 


LI2 (Lmorine diatomic) 


10.93 


2.96 


12.67 


9.55 


2.85 


13, 


.00 


10.55 


NaCl (hodium Chloride) 


8.64 


3.30 


7.55 


6.24 


3.35 


11, 


.30 


8.95 


bi(J (bilicon monoxide) 


11.60 


4.61 


12.44 


9.96 


4.74 


13, 


.48 


11.67 


Ch (Carbon monosulnde) 


11.58 


3.85 


13.59 


10.42 


3.62 


13, 


.77 


11.26 


CIO (Monocnlorine monoxide) 


8.85 


-0.03 


12.41 


8.16 


-0.26 


12, 


.58 


8.86 


Clr (Cnlorine mononuoride) 


12.43 


3.28 


14.96 


10.86 


3.08 


15, 


.57 


12.07 


■ TT fT\' '1 N 

hi2ilg (Disilane) 


11.33 


6.56 


10.09 


9.12 


6.37 


13, 


.13 


11.05 


CH3CI (Methyl chloride) 


12.01 


6.25 


12.16 


9.16 


6.28 


15, 


.05 


12.09 


CH3SH (Methanethiol) 


10.01 


4.76 


9.10 


7.53 


4.89 


11, 


.62 


10.14 


SO2 (Sulfur dioxide) 


11.74 


3.64 


13.66 


10.99 


3.51 


13, 


.95 


11.78 


Btg (Borane, tritluoro-) 


17.22 


9.55 


13.80 


12.70 


10.04 


21, 


.46 


18.10 


BCI3 (Borane, trichloro-) 


12.07 


4.85 


13.24 


10.80 


4.68 


14, 


.44 


11.31 


AICI3 (Aluminum trichloride) 


12.13 


5.66 


11.57 


9.67 


5.73 


12, 


.73 


11.27 


CI" 4 (Carbon tetraiiuoride) 


17.85 


10.68 


14.18 


13.60 


12.60 


18, 


.35 


18.12 


CCI4 (Carbon tetrachloride) 


11.97 


4.90 


11.90 


9.86 


4.86 


12, 


.23 


10.54 


OCS (Carbonyl sulfide) 


12.13 


5.51 


15.04 


10.09 


5.20 


15, 


.17 


12.69 


CS2 (Carbon disumde) 


10.19 


3.94 


12.06 


9.50 


3.74 


12, 


.83 


10.21 


Vjr 2^ (Carbonic dmuoride) 


16.08 


6.63 


18.68 


11.46 


6.39 


19, 


.13 


15.78 


Sil 4 (Silicon tetrafiuoride) 


16.95 


9.66 


14.84 


13.30 


10.87 


18, 


.38 


17.03 


N2O (Nitrous oxide) 


15.01 


6.81 


18.48 


11.89 


6.56 


18, 


.74 


15.49 


Nt 3 (Nitrogen trmuoride) 


15.76 


7.28 


21.00 


13.00 


7.17 


21, 


.65 


16.54 


PI 3 (Phosphorus trmuoride) 


13.00 


6.43 


14.13 


10.16 


6.72 


16, 


.22 


13.48 


O3 (Ozone) 


11.06 


1.72 


14.19 


10.40 


1.46 


14, 


.27 


10.95 


F2O (Difluorine monoxide) 


13.82 


3.22 


17.96 


12.34 


3.13 


18, 


.16 


13.52 


CIF3 (Chlorine trifluoride) 


11.79 


3.02 


14.48 


10.69 


3.02 


14, 


.98 


11.72 


C2F4 (Tetrafluoroethylene) 


12.45 


5.88 


16.43 


9.13 


5.72 


16, 


.80 


13.97 
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CHgGGH (Propyne) 


11.69 


6, 


.21 


9.29 


8.32 


6.49 


15, 


.16 


11.54 


CH2CCH2 (AUene) 


10.83 


5, 


.69 


13.81 


8.55 


5.55 


14, 


.75 


11.89 


C3H4 (Cyclopropene) 


11.87 


4, 


.90 


13.67 


8.13 


4.79 


14, 


.78 


11.94 


TT 1 \ 

Cgtlg (Cyclopropane) 


11.64 


6, 


.97 


9.91 


9.21 


7.64 


12, 


.86 


12.25 


CH2F2 (Methane, dmuoro-) 


14.15 


7, 


.94 


11.27 


10.37 


9.39 


18, 


.17 


14.57 


CHF3 (Methane, trifluoro-) 


15.44 


9, 


.23 


12.97 


11.94 


10.78 


17, 


.62 


16.11 


/^TT /^l /"IV /TjII 11 "iN 

CH2CI2 (Methylene chloride) 


12.18 


5, 


.86 


13.28 


9.52 


5.80 


14, 


.27 


11.91 


CHCI3 (Chloroform) 


12.38 


5, 


.31 


12.67 


9.75 


5.26 


12, 


.92 


11.20 


/^TT AT/^ /TV T J_l 'J- \ 

CH3NO2 (Methane, nitro-) 


11.94 


3 


.67 


15.16 


9.24 


3.44 


15, 


.25 


11.52 


/^TT C1*TT /TV T J_l 1*1 \ 

CHgbiHg (Methyl silane) 


12.35 


7, 


.30 


11.00 


9.87 


7.60 


13, 


.18 


12.42 


XX XX / X^ • • 1 \ 

HCOOH (Formic acid) 


11.98 


5, 


.13 


16.09 


9.03 


4.95 


16, 


.84 


13.35 


CH3CONH2 (Acetamide) 


10.05 


4, 


.98 


11.35 


7.59 


4.89 


15, 


,37 


10.91 


C2N2 (Cyanogen) 


13.90 


5, 


.83 


14.73 


12.67 


5.74 


14, 


.98 


13.28 


CH2CO (Ketene) 


10.32 


3 


.72 


14.32 


8.11 


3.42 


14, 


.70 


11.32 


C2H4O (Ethylene oxide) 


11.68 


6, 


.10 


9.16 


8.37 


7.34 


13, 


.25 


12.10 


C2H2O2 (Ethanedial) 


10.04 


1, 


.80 


11.15 


8.63 


1.57 


11, 


.37 


9.16 


CH3CH2OH (Ethanol) 


11.38 


5, 


.64 


9.37 


8.07 


6.63 


13, 


.30 


11.60 


CH3OCH3 (Dimethyl ether) 


10.79 


5, 


.47 


8.59 


7.71 


6.64 


12, 


.29 


11.31 


XX /mi \ 

C2H4S (Thiirane) 


9.93 


4, 


.34 


11.64 


7.41 


4.13 


12, 


.53 


10.12 


CH2CHF (Ethene, nuoro-) 


11.55 


5, 


.46 


14.50 


8.78 


5.41 


15, 


.41 


12.73 


CH3CH2CI (Ethyl chloride) 


11.74 


6, 


.31 


11.37 


8.90 


6.36 


15, 


.13 


11.85 


CH2CHCI (Ethene, chloro-) 


11.35 


4, 


.92 


13.61 


8.50 


4.76 


14, 


.23 


11.59 


CH3COCI (Acetyl Chloride) 


11.97 


4, 


.72 


15.12 


9.29 


4.48 


15, 


.16 


11.79 


NO2 (Nitrogen dioxide) 


9.79 


1, 


.13 


12.94 


8.85 


0.93 


13, 


.37 


10.08 


CFCI3 (Trichloromonofluoromethane) 


12.61 


5, 


.21 


12.99 


10.26 


5.11 


13, 


.32 


11.26 


r^T~^ /^i / 1\ iT j_i 11 1 • r\ \ 

CF3CI (Methane, chlorotrmuoro-) 


14.27 


7, 


.18 


16.96 


11.44 


6.89 


17, 


.67 


14.30 


HCCF (Fluoroacetylene) 


12.04 


6, 


.50 


11.49 


9.12 


6.98 


14, 


.93 


12.64 


HCCCN (Cyanoacetylene) 


12.20 


5, 


.16 


13.29 


11.03 


5.05 


13, 


.63 


12.00 


C4N2 (2-Butynedinitrile) 


11.52 


4, 


.44 


11.76 


10.28 


4.37 


11, 


.94 


10.70 


C2N2 (Cyanogen) 


13.90 


5, 


.83 


14.73 


12.67 


5.74 


14, 


.98 


13.28 
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C3O2 (Carbon suboxide) 


11.64 


4.67 


13.37 


11.42 


4.47 


13.53 


11.82 


t CM (Cyanogen fluoride) 


14.33 


7.31 


14.25 


11.41 


7.76 


17.03 


14.59 


C4H2 (Diacetylene) 


11.00 


4.77 


12.17 


9.94 


4.69 


12.67 


11.09 


HjCS (Thioformaldehyde) 


9.18 


1.80 


11.02 


8.20 


1.48 


11.27 


8.87 


CHONH2 (Formamide) 


10.81 


5.15 


15.38 


8.13 


5.00 


16.16 


11.67 


CH2CCI2 (Ethene, 1,1-dichloro-) 


11.17 


4.69 


13.08 


8.66 


4.58 


14.11 


11.01 


C2HF3 (Trifiuoroethylene) 


11.11 


5.39 


15.16 


8.62 


5.31 


15.96 


13.10 


CH2CF2 (Ethene, 1,1-difluoro-) 


11.81 


5.77 


15.06 


8.87 


5.66 


15.99 


13.12 


/^TT T~\ /"IV T J 1 in ■ 1 \ 

CH3F (Methyl fluoride) 


14.09 


7.88 


11.02 


10.21 


9.02 


17.02 


14.36 


t -1 /"T^ ' r\ I'll j1 \ 

Ci'2C'l2 (Uifluorodichloromethane) 


13.33 


6.06 


14.58 


10.86 


5.87 


15.08 


12.56 


SiF2 (Silicon difluoride) 


11.04 


4.09 


12.17 


9.62 


4.26 


13.40 


11.50 


MSE 




-7.22 


0.78 


-2.33 


-7.20 


2.73 


0.03 


MAE 




7.22 


2.11 


2.33 


7.20 


2.74 


0.47 


rms 




7.44 


2.45 


2.57 


7.48 


3.13 


0.63 



TABLE III. Comparison of errors (in eV) of the 115 derivative discontinuities (in eV) of the 



41(1 and LB94 



42| 



FG115 database \2% calculated by the DD(Oth) and DD(lst), using the LDA 
functionals. Here DD(Oth) = is obtained from the zeroth-order perturbation theory, and DD(lst), 
which is calculated by the difference between the Aks + ^xc and Aks gaps, is obtained from the 
first-order perturbation theory (i.e., "frozen orbital approximation"). The reference (DD) values 
are calculated by the difference between the Eg and Aks gaps, using the respective functionals. 
Note that since LB94 yields very accurate Eg, and possibly accurate Aks (due to the correct 
asymptote of the LB94 potential), the reference DD values calculated by LB94 should be more 
accurate than those calculated by LDA. 



Molecule 




LDA 






LB94 




DD 


DD(Oth) 


DD(lst) 


DD 


DD(Oth) 


DD(lst) 


H (Hydrogen atom) 


4.48 


-4.48 


1.86 


1.68 


-1.68 


3.10 


He (Helium atom) 


3.80 


-3.80 


0.20 


6.70 


-6.70 


0.20 


Li (Lithium atom) 


2.85 


-2.85 


1.51 


7.27 


-7.27 


1.40 


Be (Beryllium atom) 


4.01 


-4.01 


2.93 


5.58 


-5.58 


3.12 
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B (Boron atom) 


6.00 


-6.00 


4.12 


7.90 


-7.90 


3.36 


C (Carbon atom) 


8.01 


-8.01 


5.29 


10.11 


-10.11 


4.69 


N (Nitrogen atom) 


7.92 


-7.92 


2.08 


9.87 


-9.87 


3.86 


(Oxygen atom) 


9.98 


-9.98 


4.17 


12.48 


-12.48 


5.59 


F (Fluorine atom) 


12.02 


-12.02 


6.82 


14.80 


-14.80 


8.37 


Ne (Neon atom) 


6.15 


-6.15 


0.33 


8.92 


-8.92 


0.30 


Na (Sodium atom) 


2.97 


-2.97 


1.56 


6.85 


-6.85 


1.93 


TV If /TV iT * j_ \ 

Mg (Magnesmm atom) 


2.81 


-2.81 


1.95 


4.36 


-4.36 


1.67 


Al (Alummum atom) 


4.44 


-4.44 


2.20 


6.10 


-6.10 


1.58 


Si (Silicon atom) 


5.78 


-5.78 


2.61 


7.03 


-7.03 


2.40 


P (Phosphorus atom) 


5.96 


-5.96 


1.44 


7.84 


-7.84 


2.38 


S (Sulfur atom) 


7.20 


-7.20 


2.60 


8.98 


-8.98 


3.45 


CI (Chlorine atom) 


8.46 


-8.46 


4.09 


10.11 


-10.11 


3.93 


Ar (Argon atom) 


4.29 


-4.29 


0.26 


6.38 


-6.38 


0.15 


/"ITT /TV T J_l 1 1* 1\ 

CH3 (Methyl radical) 


6.24 


-6.24 


0.21 


8.67 


-8.67 


1.86 


CH4 (Methane) 


2.26 


-2.26 


0.72 


4.71 


-4.71 


0.77 


NH (Imidogen) ("'I] ) 


7.94 


-7.94 


0.68 


10.12 


-10.12 


3.60 


NH2 (Amino radical) 


7.90 


-7.90 


0.70 


10.10 


-10.10 


3.19 


NH3 (Ammonia) 


2.83 


-2.83 


0.95 


5.32 


-5.32 


1.05 


OH (Hydroxyl radical) 


9.63 


-9.63 


1.74 


12.01 


-12.01 


5.12 


TT /" \ /ttT I \ 

H2O (Water) 


3.28 


-3.28 


1.19 


6.01 


-6.01 


1.41 


(Hydrogen nuoride) 


3.50 


-3.50 


1.73 


6.70 


-6.70 


2.29 


S1H3 (Silyl) 


5.48 


-5.48 


1.68 


7.34 


-7.34 


2.42 


O'TT /O'l \ 

S1H4 (Silane) 


2.41 


-2.41 


0.81 


4.67 


-4.67 


0.58 


T~*TT /T^l 1 * \ 

PH3 (Phosphme) 


2.59 


-2.59 


0.92 


5.01 


-5.01 


2.74 


HjS (Hydrogen sulfide) 


2.88 


-2.88 


1.25 


5.50 


-5.50 


3.06 


HCl (Hydrogen sulfide) 


3.47 


-3.47 


2.51 


6.25 


-6.25 


2.59 


C2H2 (Acetylene) 


2.45 


-2.45 


5.71 


6.91 


-6.91 


2.62 


C2H4 (Ethylene) 


3.41 


-3.41 


5.41 


6.84 


-6.84 


3.12 


C2H6 (Ethane) 


2.26 


-2.26 


0.73 


4.62 


-4.62 


0.69 
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TT/^AT /TT J • J \ 

HCN (Hydrogen cyanide) 


3.66 


-3.66 


6.04 


7.46 


-7.46 


3.48 


CO (Carbon monoxide) 


6.06 


-6.06 


5.35 


8.81 


-8.81 


3.72 


HCO (Formyl radical) 


6.68 


-6.68 


4.25 


8.59 


-8.59 


3.27 


H2CO (Formaldehyde) 


5.52 


-5.52 


5.58 


8.59 


-8.59 


3.51 


CH3OH (Methyl alcohol) 


2.53 


-2.53 


1.48 


5.37 


-5.37 


1.79 


AT /AT" j_ 1 * j_ * \ 

N2 (Nitrogen diatomic) 


7.63 


-7.63 


4.67 


9.45 


-9.45 


3.64 


AT TT /TT J • \ 

M2II4 (Hydrazine) 


3.00 


-3.00 


1.05 


5.41 


-5.41 


1.06 


AT/^ /AT*J_ • • 1 \ 

NO (Nitric oxide) 


8.80 


-8.80 


5.03 


10.39 


-10.39 


4.23 


O2 (Oxygen diatomic) ( 2^^) 


9.36 


-9.36 


2.13 


10.63 


-10.63 


4.47 


TT /TT 1 • 1 \ 

H2O2 (Hydrogen peroxide) 


4.33 


-4.33 


8.71 


8.51 


-8.51 


6.39 


F2 (Fluorine diatomic) 


10.60 


-10.60 


6.94 


12.16 


-12.16 


6.11 


CO2 (Carbon dioxide) 


3.49 


-3.49 


1.53 


7.14 


-7.14 


5.33 


T~» /T^l 1 1 • j_ • \ 

P2 (r'hospnorus diatomic) 


5.66 


-5.66 


1.88 


6.50 


-6.50 


1.61 


^2 (bullur diatomic) ( Sg) 


6.35 


-6.35 


2.24 


7.20 


-7.20 


1.90 


CI2 (Chlorine diatomic) 


6.59 


-6.59 


3.12 


7.70 


-7.70 


2.45 


AT fW /C J* /~n 1 "J \ 

NaCl (bodium Chloride) 


2.94 


-2.94 


1.31 


5.60 


-5.60 


2.35 


biO (bilicon monoxide) 


5.35 


-5.35 


2.48 


6.93 


-6.93 


1.81 


Ch (Carbon monosulnde) 


6.56 


-6.56 


3.17 


7.64 


-7.64 


2.51 


/~\ /TV (T 11* * 1 \ 

UIO (Monochlorine monoxide) 


8.19 


-8.19 


4.24 


9.11 


-9.11 


3.73 


Cli' (Chlorine monoiiuoride) 


7.58 


-7.58 


4.10 


8.99 


-8.99 


3.50 


Si2Hg (Disilane) 


2.56 


-2.56 


0.98 


4.69 


-4.69 


2.08 


CH3CI (Methyl chloride) 


2.91 


-2.91 


3.01 


5.81 


-5.81 


2.96 


CH3SH (Methanethiol) 


2.76 


-2.76 


1.58 


5.25 


-5.25 


1.48 


b02 (bullur dioxide) 


7.36 


-7.36 


2.66 


8.27 


-8.27 


2.17 


T^TTI /T^ J. \ 

Bl'3 (Borane, trilluoro-) 


3.15 


-3.15 


1.10 


8.06 


-8.06 


3.35 


BCI3 (Borane, trichloro-) 


5.95 


-5.95 


2.44 


6.63 


-6.63 


3.12 


AICI3 (Aluminum trichloride) 


4.01 


-4.01 


1.89 


5.55 


-5.55 


1.46 


CF4 (Carbon tetrafluoride) 


2.91 


-2.91 


0.59 


5.52 


-5.52 


0.23 


CCI4 (Carbon tetrachloride) 


4.96 


-4.96 


2.04 


5.68 


-5.68 


1.70 


OCS (Carbonyl sulfide) 


4.59 


-4.59 


4.95 


7.49 


-7.49 


2.48 
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(Carbon disumaej 


5. 00 


-5.56 


2.5d 


0.47 


-6.47 


2.0Z 


y-jr^y^ (Carbonic dinnoridcj 


4.83 


-4.83 


7.22 


9.40 


A /I A 

-9.40 


00c 

3.35 


0*T7> /CI* J. J. £] 'J \ 

hii' 4 (Silicon tetranuoride) 


3.64 


-3.64 


1.53 


6.16 


-6.16 


1.35 


N2O (Nitrous oxide) 


5.08 


-5.08 


6.59 


8.93 


-8.93 


3.25 


Nt 3 (Nitrogen trmuoride) 


5.72 


-5.72 


8.00 


9.38 


-9.38 


5.11 


T»Tn /Til 1 J. "JH • J \ 

3 (Pnospnorus trmuoride) 


3.73 


-3.73 


3.96 


6.75 


-6.75 


2.74 


O3 (Ozone) 


8.68 


-8.68 


3.79 


9.49 


-9.49 


3.33 


(Diiiuorine monoxide) 


9.12 


-9.12 


5.63 


10.40 


-10.40 


4.64 


CIF3 (Chlorine trifluoride) 


7.67 


-7.67 


3.79 


8.70 


-8.70 


3.27 


y-j-^y ^ ( ietranuoroetnylene) 


3.25 


-3.25 


7.30 


8.25 


-8.25 


2.83 


LH3GGH (Propyne) 


2.11 


-2.11 


0.97 


5.06 


-5.06 


3.61 


CH2CCH2 (AUene) 


2.86 


-2.86 


5.26 


6.34 


-6.34 


2.86 


C3H4 (Cyclopropene) 


3.23 


-3.23 


5.55 


7.16 


-7.16 


2.83 


C3Hg (Cyclopropane) 


2.24 


-2.24 


0.70 


4.62 


-4.62 


0.61 


CH2F2 (Methane, diiiuoro-) 


2.43 


-2.43 


0.89 


5.18 


-5.18 


3.60 


CHF3 (Methane, trifiuoro-) 


2.71 


-2.71 


1.03 


5.34 


-5.34 


1.51 


CH2CI2 (Methylene chloride) 


3.67 


-3.67 


3.75 


6.11 


-6.11 


2.36 


/^TT/^l / /^l 1 f \ 

CHCI3 (Chloroform) 


4.44 


-4.44 


2.92 


5.94 


-5.94 


1.73 


CH3NO2 (Methane, nitro-) 


5.57 


-5.57 


5.92 


8.08 


-8.08 


3.73 


/"ITT O'TT /AT J.1 1 *! \ 

CH3S1H3 (Methyl silane) 


2.57 


-2.57 


1.14 


4.83 


-4.83 


0.76 


TT/^ /^TT /T~l • • 1 \ 

HCOOH (Formic acid) 


3.90 


-3.90 


7.07 


8.40 


-8.40 


3.50 


/^TT /^/^ATTT /A J- ' J \ 

CH3CONH2 (Acetamide) 


2.61 


-2.61 


3.76 


6.01 


-6.01 


4.46 


C2N2 (Cyanogen) 


6.84 


-6.84 


2.06 


7.54 


-7.54: 


1.70 


CH2CO (Ketene) 


4.39 


-4.39 


6.21 


7.90 


-7.90 


3.38 


(hjtnylene oxide) 


2.27 


-2.27 


0.79 


4.76 


-A. lb 


l.lo 


C2H2O2 (Ethanedial) 


6.83 


-6.83 


2.52 


7.59 


-7.59 


2.21 


CH3CH2OH (Ethanol) 


2.44 


-2.44 


1.29 


4.97 


-4.97 


1.70 


CH3OCH3 (Dimethyl ether) 


2.24 


-2.24 


0.88 


4.67 


-4.67 


0.99 


C2H4S (Thiirane) 


3.07 


-3.07 


4.23 


5.99 


-5.99 


2.41 


CH2CHF (Ethene, fluoro-) 


3.31 


-3.31 


5.73 


7.32 


-7.32 


2.67 
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/"ITT /^TT /^l /TT'J.I 1 11 • J \ 

CHgCHjCl (Ethyl chloride) 


2.59 


-2.59 


2.47 


5.49 


-5.49 


3.28 


CHjCHCl (Jitncnc, cmoro-j 


3.58 


-3.58 


5.10 


6.83 


-6.83 


2.64 


CH3COCI (Acetyl Chloride) 


4.57 


-4.57 


5.83 


7.31 


-7.31 


3.36 


NO2 (Nitrogen dioxide) 


7.72 


-7.72 


4.08 


9.15 


-9.15 


3.29 


CFCI3 (Trichloromonofluoromethane) 


5.05 


-5.05 


2.73 


6.15 


-6.15 


2.06 


CF3CI (Methane, chlorotrifiuoro-) 


4.27 


-4.27 


5.52 


7.41 


-7.41 


3.37 


HCCF (Fluoroacetylcnc) 


2.62 


-2.62 


2.37 


5.67 


-5.67 


2.28 


HCCCN (Cyanoacetylene) 


5.86 


-5.86 


2.26 


6.94 


-6.94 


1.63 


C4N2 (2-Butynedinitrile) 


5.84 


-5.84 


1.47 


6.34 


-6.34 


1.24 


C2N2 (Cyanogen) 


6.84 


-6.84 


2.06 


7.54 


-7.54 


1.70 


C3O2 (Carbon suboxide) 


6.75 


-6.75 


1.95 


7.35 


-7.35 


1.71 


1" CM (Cyanogen nuoride) 


4.10 


-4.10 


2.84 


6.83 


-6.83 


2.43 


C4H2 (Diacetylene) 


5.17 


-5.17 


2.23 


6.40 


-6.40 


1.58 


TT /mi ■ 111 1 \ 

H2CS (Thioformaldehyde) 


6.40 


-6.40 


2.82 


7.39 


-7.39 


2.40 


CHONH2 (Formamide) 


2.98 


-2.98 


7.25 


6.66 


-6.66 


4.49 


T T 1 / J_l 111*11 \ 

CH2CCI2 (Etliene, 1,1-dichloro-) 


3.97 


-3.97 


4.42 


6.43 


-6.43 


3.10 


TTTT^ j-1 1 \ 

C2Hi'3 (Irinuorocthylcnc) 


3.23 


-3.23 


6.54 


7.79 


-7.79 


2.86 


CH2CF2 (Ethene, 1,1-difluoro-) 


3.10 


-3.10 


6.19 


7.46 


-7.46 


2.87 


✓^TT "n /n (Til in • 1 \ 

CH3F (Methyl fluoride) 


2.33 


-2.33 


0.81 


5.34 


-5.34 


2.66 


1 — T /^l /T^ 'fl 1*11 ll \ 

Ct2Cl2 (Difluorodichloromethane) 


4.80 


-4.80 


3.72 


6.69 


-6.69 


2.52 


oir 2 (oiucon Qinuoriue ) 


0.00 


-0.00 


Z.OO 




T OA 


i.yi 


MSE 




-4.89 


3.11 




-7.23 


2.70 


MAE 




4.89 


3.11 




7.23 


2.70 


rms 




5.36 


3.73 




7.48 


3.01 
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